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Abstract 

We consider chiral U(N) x U(N) models with fermions in the limit of infinitely 
large local bare Yukawa coupling. When the scalar field is subject to non-linear con- 
straint, phase transitions in these models are seen to be identical to those in the 
corresponding purely bosonic ones. Relaxing the non-linear constraint, we compute 
the seventh-order strong-coupling series for the susceptibility in these models and 
analyze them numerically for the U(2) x U{2) case. We find that in four dimensions 
the approach to the phase transition follows to a good accuracy the mean-field criti- 
cal behavior, indicating the absence of non-trivial fixed points at strong coupling and 
being consistent with the first-order nature of the transition. In three dimensions, the 
strongly-coupled bosonic U(2) X U (2) model (without gauge fields) has a first-order 
transition strong enough to accommodate electroweak baryogenesis only for a narrow 
region of the bare parameter space. 
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Phase transitions in chiral sigma models are of interest both in three dimensions 
where they serve as models of the chiral phase transition in QCD jl], and in four 
dimensions where the nature of the phase transition determines viability of certain 
models of dynamical electroweak symmetry breaking RL HL H. In the latter case, the 
analysis should include fermions with large Yukawa couplings because they can influ- 
ence the nature of the transition. In ref. ||, it was argued, based on an approximate 
method of dealing with strongly coupled fermions, that a large Yukawa interaction 
can increase an otherwise unacceptably small hierarchy between the symmetry break- 
ing scale and the cutoff. Ref. f| considered the case of moderate (0(1)) bare Yukawa 
coupling and found no substantial increase in the hierarchy for that case. 

In this note we consider U(N) x U (N) chiral models with fermions in the limit of 
infinite bare Yukawa coupling. This limit should give us some insight into the behavior 
of systems with finite but large (much larger than one) bare Yukawa coupling y. For 
models where the order parameter is a single scalar field, it is known that the phase 
diagram in the limit y — > oo depends on the choice of a lattice form for the Yukawa 
interaction ||. This is natural when there are phase transitions that involve both 
ferromagnetic and antiferromagnetic order, because in that case the order parameter 
varies on the scale of lattice spacing and that variation cannot be removed by a 
redefinition of the field. Here we are interested in a paramagnetic-ferromagnetic 
transition, so we should choose a lattice action that has it. 

We use the local form of Yukawa interaction 

SY = yJ2(xtm + mlxt) (i) 

i 

where i labels sites of a four- dimensional lattice with unit spacing; ipf (ipf) and 
Xi {Xi) are staggered fermions, each has one spin component per site but forms a 
fundamental (conjugated fundamental) representation of the respective U (N) and in 
addition comes in N c "colors", a = 1, ...,N C . This corresponds to having 8N C Dirac 
fermions in a continuum theory. Just as in the case of a single scalar ||, in the limit 
y — > oo the fermion fields can be integrated away leaving a local correction to the 
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potential for 0. In our case, this changes scalar potential V(4>) into 

y(0) = y(0)-iV c Trln(^). (2) 

This result holds when fermions have a gauge interaction. 

In the non-linear limit of the model, when the order parameter is subject to 
the constraint <p[4>i oc 1, the correction to the potential is a constant. We see that 
bosonic correlators of a non-linear model with y — > oo coincide with those of the 
the corresponding model without fermions. Hence, the phase transitions in the two 
models are identical. If the phase transition in the model without fermions is of the 
first order, as suggested by the previous work, so is the phase transition in the model 
with fermions in the limit of infinite Yukawa coupling, and with exactly the same 
strength. (This applies also to models with an hermitean order parameter considered 
in ref.[§.) 

Even when the non-linear constraint is relaxed, there can be no second-order or 
weakly first-order ferromagnetic-paramagnetic transition in a y — > oo model if there 
is none in its bosonic counterpart. This is a simple consequence of universality, since 
the logarithmic interaction in (fj) is a local one. So, for the purpose of searching for 
non-trivial fixed points and weakly first-order transitions, it is in principle sufficient 
to study the purely bosonic model, without the logarithm, although we have also 
done independent analysis of the effective theory with the potential V. 

We have computed the strong-coupling expansion for the effective bosonic theory 
to the seventh order and done numerical analysis of the series for the U(2) x £7(2) 
case. Our results are consistent with the absence of non-trivial fixed points at strong 
coupling both in the case with and without fermions. Assuming that the phase 
transition is of the first order, we find that it is a relatively weak one for most of the 
bare parameter space. We also considered the U(2) x U (2) model without fermions 
in three dimensions and found that except for a narrow region of bare parameters 
near the line beyond which the potential becomes unbounded from below, the phase 
transition (without gauge fields) is too weakly first-order to provide deviations from 
equilibrium required for electroweak baryogenesis 0. 
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Thus, we consider the linear U(N) x U(N) model for which the bosonic action S 

is 

s = -Ert^+^l+E^) (3) 

(ij) i 

V^i) = VTr^ + A^Tr^J^AsTr^) 2 ; (4) 

<p is now an arbitrary complex 2x2 matrix. We consider only the case A2 > 0. The 
bare action @ is then bounded from below if iVAi + A2 > 0. The non-linear limit is 
recovered if we let /i 2 , A2 and iVAi + A2 to infinity with /i 2 /(iVAi + A2) fixed. 

The strong coupling expansion for (|2|) is the expansion in powers of the first term 
in (§)•[] Coefficients in that expansion are functions of the bare parameters /i 2 , Ai 
and A2 and can be expressed through linear combinations of products of invariant 
ordinary (non-functional) integrals such as 

h = Z- 1 jTr^^expl-Vi^, 

h = ^ 1 |Tr(0V) 2 exp[-y(0)]#, (5) 

h = Z- 1 |[Tr(0t0)] 2 ex P [-f(0)]#, 

h = Z' 1 /"Tr((/.V) 3 exp[-F(0)]# etc. 

where V is given in (H) and Z = j exp[— V((j))]d(j). 

The invariant integrals I n of as functions of the bare parameters, have the 
form 

I n = (2\ s )-^ 2 F n (X 2 /Xs, H 2 /^Xs) (6) 

where 

A s = Ai + A2 

and p is the total power of <p and <^ in the traces in the integrand: p — 1 for Ii, p = 2 

for I 2 and I 3 , p = 3 for I4 etc. So, for fixed A 2 /A s and /x 2 / y/2\ s , the strong coupling 

series is a series in powers of (3 = (2A S ) -1 / 2 . 

^In statistical mecnahics, such expansions are called high-temperature expansions We chose 
not to use that terminology to avoid confusion with field theories at finite temperature. 
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We have computed the seventh-order series for zero-momentum susceptibility 



X = £<Tr0o4> (7) 

i 

for three- and four- dimensional lattices in terms of the invariant integrals for a general 
U(N) x U(N) model with infinite bare Yukawa coupling. Calculation of the integrals 
and analysis of the series were done for the U(2)xU (2) case. The series were obtained 
by the recursion method of ref . [H| . The complete series will be presented elsewhere; 
here we limit ourselves to the results of the analysis.^ 

Let us consider first the purely bosonic model (N c = 0) on the four- dimensional 
simple hypercubic lattice. The values of bare parameters for which Shen || finds 
evidence for a first-order transition at large bare couplings correspond in our notation^] 
to A2/A s = 1 (Ai = 0) and /x 2 / a/2A s = 2.41775. With these parameters, the series for 
the susceptibility for the U(2) x U (2) model is 

x / Xo = 1 + 9.98596/3 + 90.66614/3 2 + 822.72732/3 3 + 7347.40485/? 4 

+ 65668.03299/3 5 + 583436.92981/3 6 + 5186731. 98528/3 7 + ... (8) 

where j3 = (2A S ) -1 ' 2 and xo is the zeroth-order susceptibility which has been factored 
out. 



The series were analyzed using Zinn- Justin's method [|T|, IB] . Denote the coeffi- 



* Comparison of our series to known limiting cases has found one discrepancy. The N = 2 

case with A2 = includes the 0(8) model. For that model on the fee (face-centered cubic) lattice 

we obtain the 7th order coefficient of 14492289.1770 while Table 8.5 of ref.flj lists 14490203.7349. 

Other coefficients through the 7th order for that model on the fee and bee (body-centered cubic) 

lattices agree. Our computer program generating the series has also reproduced, through the 7th 

order, the known results for models with fewer components of order parameter: the classical XY 

model |l^] and the general one-component model Jll| on various three-dimensional lattices, and the 

Ising model on simple and face-centered hypercubic four-dimensional lattices Jl4[ . This gives us 

confidence that our program generates series correctly. 

§ Our /i 2 is —to 2 of ref. H minus the number of nearest neighbors; our Ai and A2 are four times 

those of ref. [R| . 



5 



cients in (|8|) as a n , so that 

oo 

X/XO = E a nP n ■ (9) 



n=0 

In Zinn- Justin's method one forms the ratios 



\ n ^A 1 (10) 



(i 



2 

n-1 



Then, estimates of the susceptibility exponent 7 are obtained as 



7n = l + 2 + (11) 



and estimates of the critical "temperature" as 

/ \l/2 

P*n = — ^ exp 
For the series (13) we obtain 



(12) 



{ 73 ,... ;77 } = {1.0011,1.0013,0.9986,0.9989,0.9991}, (13) 
{P~3,..-P~j} = {9.072,9.158,8.940,8.952,8.892}. (14) 

We observe that the estimates for 7 are very close to 1 and are remarkably stable for 
such a short series. 

The values of couplings Ai and A2 to which a given estimate refers are determined a 
posteriori, through the critical value j3 c . Our estimates for j3 c and the critical exponent 
become less stable as we try to move in the region of larger /3 C or, equivalently, smaller 
X s . As a result, we could not probe the weak coupling region A s ^ 1. For the purely 
bosonic theory, the usual perturbation is applicable for small A, so this is not a very 
significant limitation. In the strong-coupling region X s > 1, except for a narrow region 
near the line 2Ai + A2 = 0, stable estimates were obtained and they were consistent 
with 7 = 1. We interpret these results as follows. 

At a second-order phase transition, the correlation length grows to infinity and 
susceptibility is singular, \ ~ (Pc — However, even for a first-order transition, 

as it is approached, the correlation length and the susceptibility may grow somewhat 
before the transition takes place. 7 = 1 is the mean-field value and it implies, in the 
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renormalization group language, that in our case the approach to the phase transition 
is controlled by the trivial fixed point Ai = A2 = 0. This fixed point is known to be 
infrared unstable by perturbation theory ]I[ - the signal of a first-order transition. 

A seventh-order series probes distances of about seven lattice spacings. If there 
were an infrared fixed point with 7 7^ 1 at strong coupling, we could detect it even 
with such a short series by taking the bare couplings close to the fixed point. Because 
non-trivial critical behavior was not found for any Ai, A2 with large A s , we can claim 
the absence of strongly-coupled infrared fixed points, except for the unlikely case that 
such a fixed point has 7 very close to 1. 

If we assume that the phase transition is, in fact, of the first order, then for 
the U(2) x U(2) model it is a relatively weak one for most of the bare parameter 
space. Indeed, we have seen that the susceptibility follows the mean-field, second- 
order, behavior up to rather large scales. Equivalently, the correlation length grows 
to the size of at least several lattice spacings. It is also instructive to compare the 
estimates (0) for the critical "temperature" a/2A s of what our series sees as a second- 
order transition with the value for which ref.@ finds a first-order transition. In our 
normalization, that valne is = 2-M = 8.94427, rather close to the numbers in 

([14]). The fact that mistaking the first-order transition for a second-order one makes 
only a small error in the critical "temperature" confirms that the transition is only 
weakly first-order. 

The correlation length for electroweak interactions is the scale of physics responsi- 
ble for the symmetry breaking and cannot be much above, say, 2 TeV. Therefore, an 
hierarchy between the correlation length and the cutoff by a factor of 10 or so should 
be sufficient for consistency of models of electroweak symmetry breaking in which the 
cutoff scale is not much above 20 TeV. 

Near the line 2Ai + A2 = 0, beyond which the potential becomes unbounded from 
below, the transition becomes more strongly first-order. We will describe the change 
in the behavior of the series in that region while discussing the three-dimensional 
case, for which this change may have an application to electroweak baryogenesis. 
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In three dimensions, our results indicate that for small A2 and large Ai, the ap- 
proach to the phase transition, at the length scale probed by the 7th order series, 
is controlled by the non-trivial 0(8) fixed point, but for most of the bare couplings 
plane, by the trivial point Ai = A2 = 0. Assumimg that the transition is of the first 
order, we find that it remains a weak one, except for the vicinity of 2Ai + A2 = 0. 
For the simple cubic lattice and A2/A s = 0.5, /i 2 / \/2A s = —5., corresponding to the 
values of fig.3 of ref.fnj, we obtain 



{ T 3,... )T7 } = {1.0005,1.0005,1.0009,1.0007,0.9997}, (15) 
{P~j,...p~j} = {1.023,1.030,1.011,1.008,1.020} (16) 

in good agreement with the mean field value 7 = 1 and the "temperature" of the 



first order transition y2\i°' « 1 deduced from fig.3 of ref. |T7[. Similar results were 
obtained for stronger couplings. This is in accord with the results of a numerical 
simulation of the non-linear £7(2) x £7(2) model in three dimensions which has 



found that to observe the first-order phase transition one has to go to lattices as large 
as 16 3 . 

Near the line 2Ai + A2 = 0, or X2/K = 2, the behavior of the series changes. 
Compare the following estimates for the critical exponent obtained for the simple 
cubic lattice, /1 2 / \/2\ s = —2. and different values of A 2 /A s : for A 2 /A s = 1.55, 1.6, 1.65, 
respectively, 

{73, ...,77} = {0.9996,0.9996,0.9990,0.9992,0.9991}, 

{ 73 , 77 } = {0.9992, 0.9993, 0.9988, 0.998 + 1.45 x 10~ 7 i, 3.16 + 0.886i} , 

{ 73 , 77 } = {2.91 + 0.218s, 5.09, -1.47 + IMi, -31.5 - 533i, 8.38 + 2.55i} . 

Imaginary values reflect the appearance of negative coefficients in the series which 
starts to display irregular behavior; the prescription for determining the signs of 
imaginary parts was chosen arbitrarily. We see that as we approach the line A2/A s = 2, 
the irregular behavior of the series sets off in lower order. In some systems, the 
presence of negative coefficients in a series signals an unphysical singularity close to 



(5 = 0, which can be mapped further away by an appropriate change of the expansion 



variable |16| . In our case, there is a physical reason for such behavior — the decrease 
in the correlation length as the phase transition becomes more strongly first-order. 

For the electroweak phase transition, the cutoff in the effective three-dimensional 
theory is of order of temperature. To prevent the washout of baryon asymmetry, the 
expectation value of the order parameter after the transition should be at least of 
order of temperature |7| and, at strong coupling, so should be the inverse correlation 
lengths in both phases. So, we find that for the strongly-coupled U(2) x U (2) model 
without gauge fields, the condition for electroweak baryogenesis is realized only in a 
narrow region near the line 2Ai + A2 = 0. 

The phase transition becomes more strongly first-order for larger U(N) x U(N) 
groups [|l 



191. However, from the results of ref. 18 we were unable to conclude 



whether for N = 3 it is strong enough to preserve the baryon asymmetry for general 
values of the couplings. 

Finally, let us turn to the y — > 00 limit of the four-dimensional U(2) x U (2) model 
with fermions. Though universality implies that there can be no infrared fixed points 
or weak first-order transitions other than those of the purely bosonic model, we have 
done an independent series analysis. The results are indeed similar to those in the 
bosonic case. Stable estimates of 7 were obtained for different N c for most of the 
strong-coupling region, X s > 1, and they were consistent with 7 = 1. The phase 
transition becomes more strongly first-order in the vicinity of the line 2Ai + A2 = 0. 
Similar results were obtained for the face-centered hypercubic lattice (with naive 
fermions), although Zinn- Justin's method often did not lead to stable estimates of 7 
in that case. 

To summarize, expanding in the inverse bare Yukawa constant y gives the effective 
bosonic theory (^) for the y — > 00 limit of chiral models of electroweak symmetry 
breaking that contain heavy fermions. The non-linear limit of that effective theory 
coincides with that of the model without fermions. For the linear U(2) x £7(2) model 
in four dimensions, the analysis of the seventh-order strong coupling series for the 
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susceptibility is consistent with the absence of any non-trivial fixed points at large 
bare self-coupling, both in the case without fermions and in the effective theory of the 
y — > oo limit. Assuming that the phase transition in the strongly-coupled U(2)xU (2) 
model is of the first order, it is a relatively weak one, both in four and (in the case 
without fermions) three dimensions, except for a region near the line beyond which 
the bare potential is unstable. 

We are grateful to T. Clark and S. Love for getting us interested in the subject and 
many helpful discussions. S.K. was supported in part by Alfred P. Sloan Foundation 
and in part by the U.S. Department of Energy under grant DE-AC02-76ER01428 
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